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We exploit results from the classical Stieltjes moment prob- 
lem to bring out the totality of all the information regarding 
phase insensitive nonclassicality of a state as captured by the 
photon number distribution p„. Central to our approach is 
the realization that n! p„ constitutes the sequence of moments 
of a (quasi) probability distribution, notwithstanding the fact 
that pn can by itself be regarded as a probability distribution. 
This leads to classicality restrictions on p„ that are local in n 
involving p„'s for only a small number of consecutive n's, en- 
abling a critical examination of the conjecture that oscillation 
in pn is a signature of nonclassicality. 

Nonclassical states of the radiation field continue to 
receive much attention. These are states for which the 
P-distribution (p{z) is not a true probability. Prominent 
among the quantitative characteristics of nonclassical- 
ity are squeezing and sub-Poissonian statistics. While 
these involve the lower order moments of ^{z), there have 
also emerged criteria involving the higher order moments. 
Among these we may note the higher order squeezing 
criteria of Hong and Mandel the related amplitude 
squared squeezing introduced by Hillery Q , and the gen- 
eralization of the Mandel Q-parameter achieved by Agar- 
wal and Tara [||. 

There has also emerged a qualitatively different kind of 
criterion for nonclassicality. While p„ = {n\p\n), which 
represents the probability of there being n photons in 
the state specified by the density operator p, is a smooth 
function of n for classical states like the coherent states 
and the thermal states, it is an oscillating function of n 
for nonclassical states like the squeezed states, as was ex- 
posed in the seminal work of Schleich and Wheeler 1^ on 
interference in phase space. Oscillation in p„ has since 
then been taken as a signature of nonclassicality [0,1). 
Indeed, these oscillations have come to be known as non- 
classical oscillations |^. It should, however, be noted 
that this oscillation criterion for nonclassicality, though 



insightful, has not been derived from basic principles and 
hence enjoys only the status of a conjecture. Its principal 
virtue lies in the fact that it is local in n, in contradistinc- 
tion to the criteria involving the Mandel Q-paramcter or 
its generalizations; the latter are expressed in terms of 
the moments of p„, and hence are global in n. 

For a radiation mode described by annihilation and 
creation operators a, measurements of operators which 
are functions of a^a (the so called phase insensitive opera- 
tors) do not depend on all the details of (p{z), but are fully 
determined by the angle averaged radial "marginal" dis- 



tribution P(/), derived from (p{z) by writing z 
and averaging over 9: 
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In particular we have 



de(p{I 



Pn = 



dIV{I)- 



r{i) = V{I)e-^ 



(1) 



(2) 



where n ~ 0,1,2,... The above relation is invertible. 
That is, the sequence {p„} represents V{I) faithfully. 

For a given state, it may happen that ^p{z) is not a true 
probability, but the phase averaged V{I) is a bonafide 
probability distribution. Such states (the Yurke-Stoler 
state pod is an example) are said to exhibit phase sensi- 
tive nonclassicality. On the other hand the nonclassical- 
ity of the state may be such that it survives the process of 
phase averaging involved in (^), thereby rendering V{I) 
itself a quasiprobability rather than a true probability. 
Then we talk of (the stronger) phase insensitive nonclas- 
sicality. Clearly, any state with sub-Poissonian statistics 
is nonclassical of the phase insensitive type. 

The purpose of this Letter is to exhibit the totality 
of all the information regarding nonclassicality of a state 
as captured by the photon number distribution sequence 
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{Pn} or, equivalently, by P{I)- Since we work at the level 
of {pn}, and not (p{z), only states with phase insensitive 
nonclassicality will be said to be "nonclassical" . All other 
states will be termed as "classical" , for brevity. 

The key to our approach is an appreciation of the 
fact that {pn} is essentially the moment sequence of a 
(quasi) probability distribution, notwithstanding the fact 
that Pn > and ^Pn = 1? Ou^d hence \j>n] can legit- 
imately be viewed as a probability distribution over the 
discrete variable n. This departure from tradition leads 
us to derive constraints on p„ which are local in n involv- 
ing p„'s for only a small number of consecutive n's, and 
enables us to critically examine the oscillation criterion in 
a direct manner. Necessary and sufficient conditions for 
absence of (phase insensitive) nonclassicality in a state 
are presented, not only in terms of the sequence {pn} 
but also in the (dual) traditional approach involving the 
factorial moments of {pn}- 

Local constraints on classical {pn}- — It turns out to be 
convenient to define a sequence {qn} in the place of {pn} 
through qn = n!p„, for n = 0,1, 2,... It follows from (y) 
that {qn} is simply the moment sequence of the distribu- 
tion P{I) = V{I)e-^: 

qn = / dip{i)r =< r >p . (3) 

Now suppose we are given a classical state so that 
^{1) > 0, for < / < cx), and consider the polynomial 
/(J) = /"(/ — x)^. Since /(/) is manifestly nonnegative 
for any real value of the parameter x, nonnegativity of 
V{I) implies, through (^), 

< /(/) >p = {x^r - 2.x/"+i + r+^)p 

= x'^qn - 2xqn+i + qn+2 > , (4) 
for all real x. That is, q„ , qn+2 > and 

g«g«+2 > ql+i , n = 0, i, 2, . . . (5) 

Written in terms of {p„}, the above condition reads 

,?l + 1 , o ,„> 

PnPn+2 > -7:)Pn+l ■ (6) 

These are our local conditions to be necessarily satisfied 
by the photon distribution {pn} of any classical state. 

Several interesting conclusions can be drawn from 
these conditions which are local in ri, and are saturated 
for every value of n by any Poissonian distribution. Sup- 
pose that we are given a state for which png = (and 
hence qno = 0) for some integer no > 0, and assume that 
the state is classical. The choice n = no in the local con- 
dition (|^) implies that qno+i — 0. Similarly the choice 
n-|-2 = ?io implies qno-i = 0- Continuing this process we 
find that for a classical state either p„ is nonzero for ev- 
ery values of n, or p„ = for all n > 0. In other words, a 



classical state other than the vacuum state, cannot be or- 
thogonal to any Fock state. To appreciate the significance 
of this conclusion, consider the state 

p = TVat^poa" , (7) 

where pQ is an arbitrary density operator, and iV is a nor- 
malization constant. We can call it the "photon added" 
Pq, for it includes the photon added coherent state |l^] 
and the photon added thermal state 1^,^ as special 
cases. Since p„ =< n\p\n >= for n < m, we con- 
clude that p is nonclassical for every m > 0. Thus, we 
have established the following result: all photon added 
states, pure or mixed, are nonclassical. 

There has been remarkable progress recently in quan- 
tum state reconstruction using techniques of optical ho- 
modyne tomography [^ . Thus, it is now possible to 
'map out' the Wigner distribution of a state using the 
inverse Radon transform, or reconstruct the density ma- 
trix in the Fock basis using a set of pattern functions. 
Schiller et al |Q report such a reconstructed p upto 
n = 6, with qn = 0.44, 0.07, 0.26, 0.30, 1.44, 3.60, 
28.80. The local conditions (||) are violated (for instance, 
0.07 X 0.30 < (0.26)2). Thus, the squeezed vacuum of 
Schiller et al turns out to be a nonclassical state of the 
phase insensitive type. That is, the nonclassicality of 
their state survives phase averaging, notwithstanding the 
fact that consideration of the Mandel Q-parameter will 
do no more than to simply indicate that this state is 
"strongly super-Poisonian" as noted by the authors. 

We now turn to the oscillation criterion for nonclassi- 
cality. There exist classical states for which the photon 
distribution {pn} is an oscillatory function of n. We may 
call these classical oscillations. Figure 1 shows such clas- 
sical oscillations for an incoherent mixture of suitably 
chosen coherent states: 
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FIG. 1. Showing classical oscillations in the photon distri- 
bution for an incoherent mixture of five coherent states with 
Xj = 0.25, 0.25, 0.2, 0.18, 0.12 and corresponding |aj|^ = 10, 
30, 60, 90, 130. Here p{n) and q{n) stand, respectively, for 
Pn and Qn of the text. Note that q„ exhibits no oscillation for 
this classical state 



On the other hand the photon added thermal (or co- 
herent) state has a {pn} with no oscillation; we have 
nevertheless seen that it is a nonclassical state. This, 
however, should not tempt one to simply dismiss the os- 
cillation criterion as being neither sufficient nor necessary 
for nonclassicality; for, as already noted, the virtue of this 
criterion does not reside in its exactitude, but rather in 
its distinction of being local in n. Thus, it is highly de- 
sirable to amend it suitably but without sacrificing its 
local character. This is easily achieved through our local 
conditions. To see this, note that (|^) implies that {qn} 
for a classical state cannot have a local maximum (if it 
had, the inequality will be violated by allowing n + 1 
to correspond to the local maximum), and hence cannot 
exhibit any oscillation. Thus we arrive at the desired 
modification: oscillation in {qn} is a sufficient condition 
for nonclassicality. It is {qn}, and not {pn}, that is the 
key to the correct oscillation criterion. Indeed, {pn} can 
oscillate for a classical state with amplitude limited by 
the extent permitted by the {n+ l)/{n + 2) factor in (||). 
Even period-two classical oscillations are allowed, as can 
demonstrated using classical states of the type (^). 

Finally, we apply our local conditions to a class of 
states obtained as superposition of two coherent states: 



1^- >= N[\zo > -fe" 



(9) 



where 9 is the relative phase (in the Pancharatnam ||T^] 
sense) between the two components of the superposition, 
and TV is the normalization constant. We have 
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(1 + (-l)"+icos( 



(10) 



It is clear that our local conditions (||) are violated by \ 
for all values of 6* 7^ ±7r/2; by odd values of n for — 7r/2 < 
d < tt/2, and by even values of n for the range —3tt/2 < 
9 < —tt/2. It is well known |l^ that coherent states are 
the only pure states for which the T'-distribution ip{z) is 
a true probability. Further, the Yurke-Stoler ||l^ states, 
which correspond to 9 = ±7r/2, have Poissonian {pn} and 
hence possess only phase sensitive nonclassicality. Thus, 
what is striking about the above analysis is the inference 
that for all values of 6* ±7r/2 the superposition state has 
phase insensitive nonclassicality, and that it is exposed 
by our lowest order local conditions! 

Necessary and sufficient condition for nonclassicality. — 
We showed that positivity of V{I) implies the local con- 
ditions (||) on its moment sequence {qn}- We now exploit 



results from the classical problem of moments to exhibit 
the necessary and sufficient conditions on {?„} in order 
that the associated state p is classical. 

The classical moment problem, on which there exists 
an enormous amount of literature [ pTf , consists of two 
parts: (i) to test if a given sequence of numbers qualifies 
to be the sequence of moments of some bonafide proba- 
bility distribution, and (ii) to reconstruct a probability 
distribution from its moment sequence. If the probability 
distribution is over the semi- infinite real line [0, 00), one 
calls it the Stieltjes moment problem. The Hamburger 
moment problem corresponds to the case where the prob- 
ability distribution is over the entire real line (— oo,cx)). 



Since, / — 



> 0, our problem of deriving the nec- 



essary and sufficient condition on the moment sequence 
{qn } in order that 7^(1) is a true probability distribution 
is indeed a Stieltjes moment problem. 

Solution of this classical problem is well known 
To exhibit this solution, construct from the moment se- 
quence {qn} two {N + l)-dimensional symmetric square 
matrices L^^^ and L^^^ as follows: 
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(11) 
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Theorem 1: The necessary and sufficient condition on 
the photon number distribution sequence {q„ = n\pn} of 
a state p, in order that the associated quasiprobability 
distribution V{I) is a true probability, is that the matri- 
ces iW, lW be nonnegative: 



L(^)>0, i(^)>0, TV = 0,1,2,. 



(12) 



It may be noted in passing that for the Hamburger 
moment problem on the entire real line (—00,00), the 
condition L'^^ > is both necessary and sufficient. 

It is immediate to relate our local condition to the 
above theorem. Nonnegativity of L'^^\L^-^'> demands, 
as a necessary condition, nonnegativity of the diagonal 
2x2 blocks of L^^\U^\ This is precisely what our 
local conditions (|^) arc! It is also clear why our local 
conditions (H) are not sufficient: positivity of the diagonal 
2x2 blocks of L'^^\L^^'' does not capture in its entirety 
the positivity of L^^-' and L^^-' required in (^2|). 

We can derive the next higher level of local conditions 
for classicality using our necessary and sufficient condi- 
tions (|l3). Given the sequence {qn}, we define 



Xn = qnqn+2/q. 



2 

n+1 ■ 



n = 0,l,2,. 



(13) 
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Then our first order local conditions (||) involving for 
three successive values of n simply reads Xn > 1, Vn, for 
any classical state. The second order local condition to 
be presented involves qn for five successive values of n or, 
equivalently, x„ for three successive values of n. 

A necessary condition for the nonnegativity of L^^\ 
L^^") is that their diagonal 3x3 blocks (such a block 
involves g„ for five successive values of n) be nonnegative. 
After some algebra this condition can be written in terms 
of the XnS as 

(x„ - 1)(:e„+2 - 1) > (^!±1^)2^ (14) 

Xn+l 

for n ~ 0, 1, 2, . . . These are our second order local con- 
ditions on {qn} or, equivalently, on {p„}. They involve 
three successive x'„s and hence five successive p'^s. Just 
like the first order conditions, these too are only nec- 
essary conditions for classicality, and we can similarly 
derive successive higher levels of local conditions. 

To see an interesting implication of (^4|), recall that a 
Poissonian distribution {pn}, i-e. a geometric sequence 
{qn}, saturates the first order local conditions and ren- 
ders a;„ = 1 identically. We now ask whether it is possible 
to have a classical state for which g„g„+2 = qn+i for some 
values of n, whereas qnqn+2 > In+i fo'^ other values of 
n. Such classical states, if they exist, can be said to be 
locally Poissonian at these former values of n. 

Suppose a classical state is locally Poissonian at some 
n = no- That is, Xn,, ~ 1. Then two applications of (14), 
once with no — n and then with no = n + 2, shows 
that the state will cease to be classical unless Xng+i — 1 
and Xno-i = 1- Continuing this process we find that 
Xn = I for all n. Thus, there exists no non-Poissonian 
classical state which is locally Poissonian: A classical 
state is either everywhere locally Poissonian (a;„ — 1 for 
all n) or is everywhere locally super-Poissonian (x„ > 1 
for all n) . 

In the light of this result we can now strengthen and 
refine our first order condition (^ by adding that for a 
classical state these inequalities are either saturated for 
all n, or they are strict inequalities for all n. 
Factorial moments. — We now present a dual approach to 
nonclassicality based on the traditional normal ordered 
moments 7„ = tr (a'''"a"p) . This approach will be seen to 
be along the lines of Agarwal and Tara However, the 
conditions for nonclassicality that we present are both 
necessary and sufficient. 

Suppose we have a state p whose normal ordered mo- 
ments 7^ (ie. factorial moments ^j^, (A;!) ^{n + ky.pn+k 
of Pn) are known. Our problem is to find necessary and 
sufficient conditions on the sequence {7™} in order that 
the state p is classical. Transcribing 7„ to the represen- 
tation in terms of the T'-distribution ip{z), and writing 
z = /^/^e'^, we have 



That is {7n} is the moment sequence of V{I), in ex- 
actly the same manner in which the sequence {qn} was 
related to ■p(i). And the state p being classical is equiva- 
lent to P{I) being a true probability distribution. Thus, 
we have a Stieltjes moment problem once again, with 
solution parallel to the earlier one. Using the moment 
sequence {7^}, form (iV -I- l)-dimensional symmetric ma- 
trices M(^),M(^) defined by 



M 



(N) 



jk 



Ij+k+l 



(16) 



where j, fc = 0, 1, . . . ,N and = 0, 1, 2, . . . 
Theorem 2: The necessary and sufficient condition that 
the state p with normal ordered (i.e. factorial) moment 
sequence {7^} be classical is that 



M(^)>0, A/(^)>0, a = 0,1, 2, 



(17) 



It should be appreciated that theorem 2 completes the 
work initiated by Agarwal and Tara |^ by improving 
their necessary condition for classicality (they had only 
the condition M^^' > 0) into the necessary and sufficient 
condition ([l7|). Thus, the constraints on the factorial 
moments {7™} arising from the requirement A/^^^ > 
are the same as in their work. However the additional 
constraints on these moments arising from the positiv- 
ity requirement on are new: with A^ = we have 

7i > 0, with A^ = 1 we have 7173 > 7I, and so on. It 
should be appreciated that these conditions cannot in- 
deed be deduced from A/^^^ > 0. 

Considering diagonal 2x2 blocks of M(^\ Af(^) we 
obtain the classicality conditions 7^7^+2 > 7^+1, for k = 
0, 1, 2, . . . Clearly, these are dual to our first order local 
conditions (|^). We may also derive conditions analogous 
to our second order local conditions ([l^), and so on. 

To conclude, either of the two approaches based re- 
spectively on the moments of P{I) and P{I) leads to 
complete solution of the problem of (phase insensitive) 
nonclassicality as coded in the photon number distribu- 
tion {pn}- It should however be appreciated that g„'s are 
well defined for every state whereas the factorial moments 
7„ may not be finite, for p„ may not decay fast enough 
as a function of n. For the states for which 7„ exists for 
all n, the two approaches are equivalent. Even then, it is 
unlikely that the connection between nonclassicality and 
oscillation in p„ could have been so easily settled in terms 
of 7„. Finally, the first approach in terms of local condi- 
tions of Pn has a distinct advantage at least in situations 
where, for some reason, p„ is known not for all values of 
n. The density matrix from Ref. |lj] which we analysed 
is such an example. The point is, even with knowledge 
of Pn only for a finite set of values of n one can now look 
for signatures of nonclassicality. 
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dip{i)r = . 



(15) 
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